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Abstract 


Measurements of diffractive dissociation cross sections in pp collisions at = 7 TeV 
are presented in kinematic regions defined by the masses Mx and My of the two 
final-state hadronic systems separated by the largest rapidity gap in the event. Differ¬ 
ential cross sections are measured as a function of = ^x /s in the region —5.5 < 
logjQ ^x < —2.5, for log^g My < 0.5, dominated by single dissociation (SD), and 0.5 < 
logjQ My < 1.1, dominated by double dissociation (DD), where Mx and My are given 
in GeV. The inclusive pp cross section is also measured as a function of the width of 
the central pseudorapidity gap Arj for Arj > 3, log^g Mx > 1.1, and log^g My > 1.1, a 
region dominated by DD. The cross sections integrated over these regions are found 
to be, respectively, 2.99 ± 0.02 (stat)lQ 29 (syst)mb, 1.18 ± 0.02 (stat) ± 0.13 (syst) mb, 
and 0.58 ± 0.01 (stat)]'^Q (syst) mb, and are used to extract extrapolated total SD and 
DD cross sections. In addition, the inclusive differential cross section, da/dAt]^, for 
events with a pseudorapidity gap adjacent to the edge of the detector, is measured 
over Arj^ = 8.4 units of pseudorapidity. The results are compared to those of other 
experiments and to theoretical predictions, and found compatible with slowly-rising 
diffractive cross sections as a function of center-of-mass energy. 
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1 Introduction 

A significant fraction (%25%) of the total inelastic proton-proton cross section at high ener¬ 
gies can be attributed to diffractive interactions, characterized by the presence of at least one 
non-exponentially suppressed large rapidity gap (LRG), i.e. a region of pseudorapidity rj de¬ 
void of particles, where for a particle moving at a polar angle 9 with respect to the beam 
J] = — ln[tan(0/2)]. If this rj region is adjacent to the diffractively scattered proton it is called 
a forward pseudorapidity gap. In hadronic interactions an LRG is presumed to be mediated 
by a color-singlet exchange carrying the vacuum quantum numbers, commonly referred to as 
Pomeron exchange. Figurej^defmes the main types of diffractive processes: single dissociation 
(SD), double dissociation (DD), and central diffraction (GD). 
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Figure 1: Schematic diagrams of (a) nondiffractive, pp —> X, and diffractive processes with (b) 
single dissociation, pp ^ Xp or pp —> pY, (c) double dissociation, pp —> XY, and (d) central 
diffraction, pp —> pXp; X(Y) represents a dissociated proton or a centrally produced hadronic 
system. 


Inclusive diffractive cross sections cannot be calculated within perturbative quantum chromo¬ 
dynamics, and are commonly described by models based on Regge theory (see e.g. Ill and ref¬ 
erences therein). The predictions of these models generally differ when extrapolated from the 
Tevatron center-of-mass energies of < 1.96 TeV to LHG energies. Therefore, measurements 
of diffractive cross sections at 7 TeV provide a valuable input for understanding diffraction and 
improving its theoretical description. They are also crucial for the proper modeling of the full 
final state of hadronic interactions in event generators, and can help to improve the simulation 
of the underlying event, as well as of the total inelastic cross section. 

The DD cross section has been recently measured at v^s = 7 TeV by the TOTEM collabora¬ 
tion |2|, for events in which both dissociated-proton masses are below ~12GeV. Other mea¬ 
surements of diffractive cross sections at the LHG, with higher dissociation masses, have ei¬ 
ther a limited precision |3l or no separation between SD and DD events lH. In this paper, we 
present the first GMS measurement of inclusive diffractive cross sections at \/s = 7 TeV. This 
measurement is based on the presence of a forward LRG, with SD- and DD-dominated event 
samples separated by using the GASTOR calorimeter 15], covering the very forward region, 
—6.6 < Tj < —5.2. A data sample with a central LRG, in which DD dominates, is also used. 
In addition, the inclusive differential cross section, dtr/ dArj^, for events with a pseudorapidity 
gap adjacent to the edge of the detector, is measured over Arj^ = 8.4 units of pseudorapidity, 
and compared to a similar ATLAS measurement |4]. The results presented here are based on 
the first GMS data collected at \/s = 7 TeV during the 2010 LHG commissioning period, when 
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3 Monte Carlo simulation 


the probability of overlapping pp interactions in the same bunch crossing (pileup), which may 
spoil the detection of the gap, was low. 


The paper is organized in eleven sections and two appendices. The CMS detector is described 
in Section Section presents the Monte Carlo (MC) simulations used in the analysis. The 
event selection and the diffractive event topologies used to measure the cross sections are dis¬ 
cussed in Sections and respectively. Sections and present the measurement of the 
forward-gap and central-gap differential cross sections, and the integrated cross sections, re¬ 
spectively. All cross sections are extracted within the detector acceptance and with minimal 
model-dependent systematic uncertainties. The extrapolation of the measured cross sections 
to the low mass regions is discussed in Section]^ Section [^presents the measurement of the 
pseudorapidity gap cross section and its comparison to the ATLAS result (H. The systematic 
uncertainties for all the measurements are discussed in Section A summary is given in 


Section 12 Appendices A and B present additional comparisons between the diffractive MC 
models used. 


2 The CMS detector 

A detailed description of the CMS detector, together with a definition of the coordinate system 
used and the relevant kinematic variables, can be found in Ref. |0. The central feature of the 
apparatus is a superconducting solenoid of 6 m internal diameter, providing a 3.8 T axial field. 
Within the field volume are located a silicon pixel and strip tracker, a crystal electromagnetic 
calorimeter (ECAL), and a brass/scintillator hadron calorimeter. Muons are measured in gas- 
ionization detectors embedded in the steel flux-return yoke of the magnet. The calorimeter 
cells are grouped in projective towers, of granularity Arj x Acp = 0.087 x 0.087 at central ra¬ 
pidities and 0.175 x 0.175 at forward rapidities. In addition to the barrel and endcap detectors, 
CMS has extensive forward calorimetry. The forward component of the hadron calorimeter, 
HE (2.9 < \ri\ < 5.2), consists of steel absorbers with embedded radiation-hard quartz fibers, 
providing fast collection of Cherenkov light. The very forward angles are covered at one end 
of CMS (—6.6 < rj < —5.2) by the CASTOR calorimeter [|5]|, made of quartz plates embedded 
in tungsten absorbers, segmented in 16 (j>-sectors and 14 z-modules. 

Two elements of the CMS monitoring system, the Beam Scintillator Counters (BSC) and the 
Beam Pick-up Timing experiment (BPTX) devices, are used to trigger the CMS readout. The 
two BSC are located at a distance of ±10.86 m from the nominal interaction point (IP) and are 
sensitive in the \ tj\ range 3.23 to 4.65. Each BSC consists of 16 scintillator tiles. The BSC elements 
have a time resolution of 3 ns and an average minimum ionizing particle detection efficiency of 
96.3%. The two BPTX devices, located around the beam pipe at a distance of 175 m from the IP 
on either side, are designed to provide precise information on the bunch structure and timing 
of the incoming beams, with better than 0.2 ns time resolution. 


3 Monte Carlo simulation 

Monte Carlo simulations are used to correct the measured distributions for the geometrical 
acceptance and reconstruction efficiency of the CMS detector, as well as for migrations from 
true to reconstructed values in the distributions of kinematic variables. We use PYTHIA 8.165 m 
|8l to generate samples of inelastic events. We compare the detector-level data distributions to 
the PYTHIA 8 4C IjHl and PYTHIA 8 MBR (Minimum Bias Rockefeller) |9| simulations and extract 
integrated cross sections using PYTHIA 8 MBR. 
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Diffractive events in the pythia 8 4C simulation are generated according to the Schuler-Sjdstrand 
model implemented in PYTHIA 6 fT], The 4C tune |8l includes a downward scaling of the 
Schuler-Sjostrand SD and DD cross sections at \/s = 7TeV by 10 and 12%, respectively 

The PYTHIA 8 MBR generator predicts the energy dependence of the total, elastic, and inelastic 
pp cross sections, and fully simulates the main diffractive components of the inelastic cross 
section: SD, DD and CD. The diffractive-event generation in pythia 8 MBR is based on a phe¬ 
nomenological renormalized-Regge-theory model UTOUTTlI , which is unitarized by interpreting 
the Pomeron flux as the probability for forming a diffractive rapidity gap. The model was orig¬ 
inally developed for the CDF experiment at the Tevatron and has been successfully tested with 
the CDF results on diffraction. The pythia 8 MBR simulation assumes a linear parametriza- 
tion of the Pomeron trajectory, a(f) = 1 -|- e -|- a'f, where t is the square of the four-momentum 
transfer between the two incident protons. We use a' = 0.25 GeV^^, and e = 0.08 IfT^ or 
£ = 0.104 lll4llT^ , to account for the possible energy dependence of e in the range of diffrac¬ 
tive masses accessible in this analysis. We find that the simulation with e = 0.08 gives a good 
description of the data. Scaling the DD cross section downwards by 15%, which preserves the 
agreement with the CDF data, further improves the description of the DD-dominated data at 
\/s = 7TeV. These modifications are incorporated into the simulation used here. The mea¬ 
sured cross sections are also compared to the predictions of pythia 6 Z2* Ifl^ and to MC gen¬ 
erators based on Regge-Gribov phenomenology: PHOJET (version 1.12-35) IITTIITSI . QGSJET-II 
(versions 03 and 04) IIT^l20ll . and EPOS LHG (HI; the latter two are commonly used in cosmic- 
ray physics (22] . 

At the stable-particle level (where stable particles are those with lifetime t such that ct > 

10 mm), the kinematic regions covered by the present measurements are defined by the masses 
Mx and My of the two final-state hadronic systems separated by the largest rapidity gap in 
the event. For a final-state particle of energy E and longitudinal momentum pz, rapidity is 
defined as y = (1/2) ln[(£ -|- pz)/(E — pz)]. At stable-particle level the gap is defined as the 
largest rapidity separation between stable particles, without any acceptance restriction. The 
final state is then separated into systems X and Y, which populate the regions on the positive 
and negative side of the rapidity gap, respectively. The corresponding masses Mx and My are 
calculated from the full set of four-vectors in the respective group of stable particles. In the 
following Sections, Mx and My are given in units of GeV. 

We use the pseudorapidity variable to select diffractive events at the detector level. At the 
stable-particle level, the true rapidity is used. For the pseudorapidity gap cross section, pseu¬ 
dorapidity (not true rapidity) is used at the hadron level to avoid unnecessary large bin mi¬ 
grations between the distributions measured at the detector and stable-particle levels. As the 
central GMS detector is insensitive to low-mass diffraction, we use the pythia 8 MBR simula¬ 
tion, which describes the data well, to extrapolate the measured cross section into the low mass 
region. 

The detailed MG simulation of the GMS detector response is based on Geant 4 (23l . Simulated 
PYTHIA 8 4G and pythia 8 MBR events are processed and reconstructed in the same manner 
as collision data. 


4 Event selection 

The present analysis is based on event samples collected during the 2010 commissioning pe¬ 
riod, when the LHG was operating at low pileup. For the results presented in Sections |5j|^ 
only data with information from the GASTOR calorimeter are used, which correspond to an 
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4 Event selection 


integrated luminosity of 16.2 and have an average number of inelastic pp collisions per 
bunch crossing of /r = 0.14. The results based on pseudorapidity-gap events presented in Sec¬ 
tion 10 are extracted from a different set of data taking runs with negligible pileup = 0.006) 

-1 


that correspond to an integrated luminosity of 20.3 Events were selected online by requir¬ 
ing a signal in both BPTX detectors, in conjunction with a signal in any of the BSC scintillators. 
These conditions correspond to requiring the presence of two crossing bunches along with ac¬ 
tivity in the main CMS detector (minimum bias trigger). 
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Figure 2: Detector-level distributions of the energy of PF objects in four pseudorapidity inter¬ 
vals: \i]pf\ < 1-4, 1.4 < \i]pf\ < 2.6, 2.6 < \t]pF\ < 3.2, and \^pf\ > 3.2, corresponding to 
the barrel, endcap, endcap-forward transition, and forward detector regions (columns), for five 
particle candidate types: charged hadrons (tracks), photons, neutral hadrons, and two types 
that yield electromagnetic or hadronic energy deposits in HF (rows). Electron and muon can¬ 
didates constitute less than 0.1% of the PF objects reconstructed in the \i^pf\ < 2.6 region, and 
are not shown. The data are compared to the predictions of the PYTHIA 8 MBR simulation, 
normalized to the integrated luminosity of the data sample. The contribution of each of the 
generated processes is shown separately. 


Offline selections Il24l are applied to remove beam-scraping, beam-halo, and noise events. In 
addition, a minimal activity in the main CMS detectors is imposed offline by requiring at least 
two particle-flow (PF) objects reconstructed within the geometrical acceptance of the BSC de¬ 
tectors (3.23 < \r]\ < 4.65, with an energy of at least 4GeV for each PF object). Particle-flow 
objects Il25ll26l are particle candidates obtained by optimally combining the information from 
the tracking system and the calorimeters. In the forward regions (|?/| > 2.5), where there is no 
tracking, PF objects reduce to calorimeter towers. Figure|^shows the distributions of the energy 
of the PF objects reconstructed in different detector regions for different particle candidates, 
compared to the prediction of the PYTHIA 8 MBR simulation, which describes the data well. 
The requirement on the minimum energy of PF objects was found by studying data collected 
in dedicated runs with no beam; it depends on the detector region and the particle candidates, 
and varies from zero for tracks to 4 GeV for the HF towers. To assure a reliable Monte Garlo 
description of the data, the two innermost (most forward) rings of HF are not used in the anal- 
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Figure 3: Event topologies in final-state particle rj space. Detector level: nondiffractive events 
(ND), diffractive events with a forward pseudorapidity gap on the positive (FGl) or negative 
(FG2) Tj-side of the detector, or with a central pseudorapidity gap (GG). Generator level: (a) 
ND, pp —> X, (b) SDl, pp ^ Xp, (d) SD2, pp —> pY, and (c, e ,f) DD, pp —> XY, events. The 
empty box represents the central GMS detector i\f]\ ^ 4.7), filled full boxes indicate final-state 
hadronic systems or a proton—the vertical thin bar at the right/left end of sketch (b)/(d). The 
dotted empty boxes in (d) and (e) represent the GASTOR calorimeter (—6.6 < rj < —5.2). 


ysis, thus limiting the central GMS detector coverage to \rj\ < 4.7. The two outermost (most 
central) HF rings are also not used for the same reason. No vertex requirement is imposed. 
This procedure gives high acceptance for diffractive events with the hadronic system outside 
the tracking acceptance (i.e. with low to moderate diffractive masses, 12 < Mx < 100 GeV). 
According to the PYTHIA 8 MBR simulation, the selection described above accepts about 90% 
of the events corresponding to the total inelastic cross section in the region of log^g Mx >1.1 
or log^Q My > 1.1. 


5 Diffractive event topologies 

The events satisfying the selection described in Section constitute a minimum bias sample 
dominated by inclusive inelastic events in the region covered by the central GMS detector (| ?/1 < 
4.7). They are mostly composed of nondiffractive (ND) events for which final-state particle 
production occurs in the entire t] space available, as shown schematically in Fig. |^. In contrast, 
diffractive events are expected to have an LRG in the final state. Experimentally, the following 
diffractive topologies are defined, depending on the position of the reconstructed LRG in the 
central detector: 
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5 Diffractive event topologies 
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Figure 4: Detector-level distributions for the (a) ?/max/ (b) ^min, and (c) Arj^ = ~ Vmin vari¬ 

ables measured in the minimum bias sample (with only statistical errors shown), compared to 
predictions of the PYTHIA 8 MBR simulation normalized to the integrated luminosity of the 
data sample. Contributions from each of the MC-generated processes, and simulated events 
with at least two overlapping interactions of any type (pileup), are shown separately The 
dashed vertical lines indicate the boundaries for the ffmax < 1/ ^mm > — 1/ and Arj^ > 3 selec¬ 
tions. 
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• FGl: a forward pseudorapidity gap at the edge of the detector on the positive rj-side 
(Figs.|^,c); 

• FG2: a forward pseudorapidity gap at the edge of the detector on the negative rj-side 
(Figs.|^,e); 

• GG: a central pseudorapidity gap in the detector around t] = 0 (Fig.|^). 

The experimental topology with a forward pseudorapidity gap on each edge of the detector 
(GD topology, Fig.j^) is neglected in this analysis because of the limited number of such events. 
For the FGl and FG2 topologies the pseudorapidity gap is related to the variables ?/max and ?/min 
(Fig.0, defined as the highest (lowest) rj of the PF object reconstructed in the central detector. 
Experimentally, the pseudorapidity gap in GG events may be expressed as Arj^ = ~ Vmin' 

where iVmin) closest-to-zero rj value of the PF objects reconstructed on the positive 
(negative) tj-side of the central detector (Fig.|^. 

Figure 1^ shows the distributions of rj^^ax, ^min, and “ Vmin minimum bias 

sample defined in Section|^ compared to MG predictions. For the Arj^ selection, the additional 
requirement that activity be present on both rj-sides of the central detector is imposed. The 
data are dominated by the contribution from ND events, for which rapidity gaps are exponen¬ 
tially suppressed Il27l . Diffractive events appear as a flattening of the exponential distributions, 
and dominate the regions of low J/max, high J/min/ and high Arj^. The absence of events around 
l^maxKl^minl) ^ 3 in Fig.|^ l^) reflects the fact that the two outermost (most central) rings of 
HF are not used in the analysis; the depletion of events around | rjmax \ (| Vmin |) ~ 2.4 corresponds 
to the transition region between the tracker and the forward calorimeters, where higher thresh¬ 
olds are applied for the latter. The regions of 3 < A//° < 6 and A?/° ^ 6 in Fig. correspond 
to the DD topology for which one or both of the ^nd edges are in the HF calorime¬ 
ters. In order to select samples of FGl, FG2, and GG events with a central LRG signature, the 
requirements rjmax < 1/ ^min > — 1/ and Arj^ > 3 are imposed, respectively. 

According to the expectations of the PYTHIA 8 MBR simulation, the event samples defined ex¬ 
perimentally as FGl or FG2 (Figs. |^,b) originate from approximately equal numbers of SD 
events with 1.1 < log^gMx ^ 2.5 and DD events for which one dissociated-proton mass is 
in this Mx range, while the other is small and escapes detection in the central detector, cf. 
Figs. |^,e. For the FG2 topology, GASTOR (—6.6 < rj < —5.2) is used to separate diffractive 
events into two samples: log^gMy < 0.5 (SD enhanced) and 0.5 < log^gMy < 1.1 (DD en¬ 
hanced, Figs.|^,e). The detection of the low-mass dissociated system, Y, is performed by using 
a GASTOR tag, defined as the presence of a signal above the energy threshold (1.48 GeV) in at 
least one of the 16 (^-sectors of the first five GASTOR modules. Since no detector is available 
for tagging the low-mass dissociated system on the positive tj-side, the FGl sample is treated 
as a control sample in this analysis. 

The range of the dissociation mass Mx for the true SD process in the FG2-type sample af¬ 
ter all detector selections is shown as a hatched histogram in Fig. for PYTHIA 8 MBR (left) 
and PYTHIA 8 4G (right), and corresponds to 1.1 < log^gMx < 2.5. Similar distributions are 
obtained for events in the FGl-type sample, in which the dissociated system originates from 
the proton on the other side of the detector. The ranges of dissociation masses, Mx and My, 
for the true DD events in the minimum bias sample after the trigger selection, in the FG2- 
type sample with a GASTOR tag, and in the GG-type sample after all detector selections, are 
shown in the efficiency plots of Figs. |^, |^, and |^, respectively. The FG2-type events, with 
the pseudorapidity gap reconstructed at the edge of the central detector, populate the region 
of 1.1 < log^gMx < 2.5 and 0.5 < log^^gMy < 1.1 (solid box in Fig. |^). The selection based 
on the GG topology requires both diffractive masses to be in the central detector; this leads to 
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6 Forward pseudorapidity gap cross sections from the FG2 event sampie 


different coverage in the (Mx, My) plane. Events populate the region of log^gMx ^ 1-1 and 
log^Q My > 1.1 (Fig.|^), in addition to Atj^ > 3, thus providing a complementary measurement 
of the DD cross section. 
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Figure 5: Simulated distributions of the dissociated mass Mx at stable-particle level for the 
SD process in the FG2 sample at successive selection stages (trigger, minimal detector activity 
within BSC acceptance, J/min > —1) for pythia 8 MBR (left) and pythia 8 4C (right). The MC 
samples are normalized to the luminosity of the data sample. 
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Figure 6: Simulated (pythia 8 MBR) event selection efficiency in the Mx vs. My plane for true 
DD events after (a) the trigger selection, and (b) the FG2 selection with a CASTOR tag or (c) the 
CG selection (Fig. |^. The regions delimited by the solid (red) lines in (b) and (c) are those of 
the cross section measurements; the dashed (red) box in (b) corresponds to the enlarged region 
for which the cross section is given (Section]^, assuming the same dependence on Mx and My; 
the dashed (blue) line in (c) marks the region of Arj > 3. 


6 Forward pseudorapidity gap cross sections from the FG2 event 
sample 

The forward pseudorapidity gap cross sections are measured as a function of the variable 
which is related to the mass Mx of the dissociated system by: 


= 



( 1 ) 


For the FG2 sample, Mx corresponds to the dissociated system that can be detected in the 
central detector (right-hand side of Figs.|^,e). The CASTOR calorimeter allows the detection 
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of the hadronic system Y when it escapes the central detector, and the separation of the FG2 
sample into subsamples corresponding to log^gMy < 0.5 and 0.5 < log^gMy < 1.1, which 
are dominated by SD and DD events, respectively. For the purely SD events, represents the 
fractional longitudinal momentum loss of the incoming proton. 

At detector level, the variable is reconstructed as 



LiiPTpi) 

V~s 


( 2 ) 


where i runs over all PF objects measured in the central detector, and and ^re the energy 
and the longitudinal momentum of the zth PF object, respectively. The energy is related to the 
particle three-momentum assuming a mass that depends on the PF object type; e.g. for charged 
hadrons a pion mass is assumed. The signs (±) in Eq. (|^ indicate whether the dissociated 
system is on the ±z side of the detector. For the FG2-type events under study, ^x corresponds 
to^+. 
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Figure 7: Two-dimensional distribution of reconstructed vs. generated ^x values for the 
events in the SD2 sample obtained with the PYTHIA 8 MBR simulation. The solid red line 
represents the condition log^Q = log^^Q ^x- 

Since part of the hadronic system X escapes the detector through the forward beam hole, and 
since low-energy particles remain undetected because of the PF object thresholds, the recon¬ 
structed values are underestimated. This can be seen in Fig. which shows a scatter plot 
of reconstructed vs. generated values of ^x for PYTHIA 8 MBR events in the FG2 sample. As 
^x decreases, its unmeasured fraction increases (the beam hole size is fixed), resulting in a 
larger deviation from the log^^Q = log^Q ^x line. The calibration factor C(^^), which brings 
the reconstructed values of (Eq. (j^) to their true values (Eq. Q) according to the formula 
log^g = log^Q -|- C(^^), is evaluated from the PYTHIA 8 MBR simulation, by studying 
the log^g ^x ~ log^Q difference in bins of log^g The factor C(^^) decreases from the value 
of 1.1 at log^g ~ —6.5 to 0.2 at log^g ~ —2.5, with an uncertainty of 9%, estimated by 
comparing the PYTHIA 8 MBR and PYTHIA 8 4G simulations. 

Figure [^presents the distribution of the calibrated log^g for the FG2 sample, compared to 
the predictions of the PYTHIA 8 MBR (top) and PYTHIA 8 4G (bottom) simulations. Figure 
shows the comparison for the entire FG2 sample. The separation of the SD and DD processes 
(hatched green and solid yellow histograms, respectively) by means of the GASTOR tag is 
clearly seen in Figs.|^,c. Overall, the PYTHIA 8 MBR MG describes the data better than PYTHIA 
8 4G and is therefore used to extract the diffractive cross sections. Both MG predictions are 
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6 Forward pseudorapidity gap cross sections from the FG2 event sampie 


presented for a Pomeron trajectory with e = 0.08 and describe the region of low well. At 
higher values, e = 0.104 would be more appropriate Ifl4l . providing a better agreement 
with the data in that region. Since the Mx dependence of e is currently not available in the MC 
models used for this analysis, we extract cross sections using e = 0.08 and evaluate a systematic 
uncertainty related to the possible variation of e, as explained in Section]^ 


16.2 Mb-’(7 TeV) 16.2 Mb ' (7 TeV) 16.2 Mb ’(7 TeV) 



Figure 8: Detector-level distributions of the reconstructed and calibrated ^x for (a) the entire 
FG2 sample, and the FG2 subsamples with (b) no GASTOR tag, and (c) a GASTOR tag (statis¬ 
tical errors only). The data are compared to the predictions of the PYTHIA 8 MBR (top three 
plots) and PYTHIA 8 4G (bottom three plots) simulations, which are normalized to the inte¬ 
grated luminosity of the data sample. The contribution of each of the generated processes is 
shown separately. 


The differential cross sections measured in bins of ^x, separately for log^Q My < 0.5 (no GAS¬ 
TOR tag) and 0.5 < log^Q My <1.1 (GASTOR tag), are calculated with the formula 


da _ 

dlogjo^x i:(Alogio^x)bin' 


( 3 ) 


where is the number of events in the bin, corrected for acceptance and migration effects, C 
is the integrated luminosity, and (A log^g ^x)bin is the bin width. The acceptance and migration 
corrections are evaluated with the iterative Bayesian unfolding technique I28i , as implemented 
in the ROOUnfOLD package [29l, with four iterations. The number of iterations is optimized 
following the procedure suggested in Ref. Il28l . by studying the difference in (goodness-of- 
fit) values after consecutive iterations; the final unfolded distribution, folded back to the detec¬ 
tor level, is consistent with the observed data. The response matrix is obtained with PYTHIA 
8 MBR (e = 0.08). The results are presented in Table they include a small correction for 
overlapping pp collisions (~7%), evaluated by comparing MG simulations with and without 
pileup. 

Figures |^,c present the measured cross sections compared to the PYTHIA 8 MBR, PYTHIA 8 4G, 
and PYTHIA 6 Z2* predictions. The error bars of the data points are dominated by systematic 
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Table 1: The differential forward pseudorapidity gap cross sections dtr/dlog^Q^x for 
log^gMy < 0.5 (SD dominated, without CASTOR tag) and 0.5 < log^gMy < 1.1 (DD dom¬ 
inated, with CASTOR tag). The first and the second errors correspond to statistical and sys¬ 
tematic uncertainties, respectively. 


bin dCTno-CASTOR/d log (mb) dcrcASTOR/d log ^x (mb) 


-5.5 

< 

logio 


< 

-5.0 

1.17 

± 

0 02+°’^® 
o-o^_0.ii 

0.30 ±0.01^“;“^ 

-5.0 

< 

logio 


< 

-4.5 

1.16 

± 

o-o^_0.i7 

0.30 ± 0.01 ± 0.04 

-4.5 

< 

logio 

?x 

< 

-4.0 

0.91 

± 

o-o^_0.l2 

0.26 ± 0.01 ± 0.03 

-4.0 

< 

logio 

?x 

< 

-3.5 

0.88 

± 

o-o^_0.09 

0.32±0.0llg;g| 

-3.5 

< 

logio 

?x 

< 

-3.0 

0.98 

± 

O-O^_0.13 

0.51 ±0.01^^;^^ 

-3.0 

< 

logio 

?x 

< 

-2.5 

0.78 

± 

o-oj-o.og 

0.67±0.03^°;J2 


uncertainties, which are discussed in Section 11 The predictions of PYTHIA 8 MBR are shown 
for two values of the e parameter of the Pomeron trajectory. Both values, e = 0.08 and e = 0.104, 
describe the measured cross section for log^g My < 0.5. The data for 0.5 < logjg My < 1.1 favor 
the smaller value of £, specifically in the region of lower ^x, corresponding to the topology in 
which both dissociation masses are low. The prediction of the Schuler-Sjostrand model, used in 
the PYTHIA 8 4C simulation, describes well the measured cross section for 0.5 < log^g My < 1.1, 
while the PYTHIA 6 Z2* simulation overestimates it. Both predictions are higher than the data 
for log^Q My < 0.5 at high log^g ^x, and the predicted rising behavior of the cross section is not 
confirmed by the data in the region of the measurement, —5.5 < log^g ^x < —2.5. 


Figures |^,d present a comparison of the measured cross sections with the PHOJET, QGSJET- 
II 03, QGSJET-ll 04, and EPOS predictions. None of the models is able to describe the magnitude 
of the cross section in the region 0.5 < log^gMy < 1.1. For log^gMy < 0.5, the PHOJET and 
EPOS generators fail to describe the falling behavior of the data, QGSJET-ll 03 describes the 
measured cross section reasonably well, while QGSJET-ll 04 underestimates the magnitude of 
the cross section. 


7 Central pseudorapidity gap cross section from the CG event 
sample 

The cross section for events with a central pseudorapidity gap is measured as a function of the 
variable A?/, defined as At] = — log^, where ^ = M^-My/(swip), with nip the proton mass. 
For purely DD events, the position of the gap center is related to the dissociation masses by the 
expression rjc = log(My/Mx). 

As discussed in Section]^ the central-gap width (Fig.|^) is reconstructed as Arj^ = ~ Vmin- 

The calibration factor C, which corrects Arj^ for detector effects according to the formula = 

A/jjgc “ C, is extracted from the PYTHIA 8 MBR MC as the difference C = Ajj^gg — Ai/gg^,. It 
amounts to C = 2.42 ± 0.12, with the uncertainty estimated from a comparison with PYTHIA 

8 4C. Figure presents the distribution of the calibrated Arj^ for the CG sample along with 
simulated distributions from PYTHIA 8 MBR and PYTHIA 8 4C. 

The differential cross section, measured in bins of A?/for A?/ > 3, logjgMx > l.l,andlogjQMy > 
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7 Central pseudorapidity gap cross section from the CG event sample 



Figure 9: Cross sections dcr/d log^Q for (a,b) log^g My < 0.5 (SD dominated) and (c,d) 
0.5 < log^gMy < 1.1 (DD dominated) compared to MC predictions: (a,c) PYTHIA 8 MBR, 
PYTHIA 8 4C, PYTHIA 6 Z2*, and (b,d) PHOJET, QGSJET-II 03, QGSJET-II 04, EPOS. Error bars are 
dominated by systematic uncertainties (discussed in Section [TT||. 


1.1, is calculated according to the formula 


dcr _ 

dAtj /3(A/y)bjn 


(4) 


where N®'"* is the number of events in a given bin, corrected for acceptance and migration 
effects, and also for the extrapolation from At]^ > 3 (for gaps overlapping rj = 0) to Arj > 3 
for all gaps, C is the integrated luminosity, and (A)y)bin is the bin width. The acceptance and 
migration corrections are evaluated with the iterative Bayesian unfolding technique Il28l with 
two iterations, optimized as described in Section The response matrix is obtained using 
PYTHIA 8 MBR with e = 0.08. 


The measured differential cross section is presented in Table and compared to predictions 
of theoretical models in Fig. 11 The results take into account the pileup correction, and the 


uncertainties are dominantly systematic (see Section]^. The prediction of the PYTHIA 8 MBR 
MC simulation with both e = 0.08 and 0.104 describes well the central-gap data. PYTHIA 8 4C 
underestimates the data in all bins, while PYTHIA 6 Z2* overestimates the data in the lowest At] 
bin. The PHOJET, QGSJET-II 03, QGSJET-II 04, and EPOS generators underestimate the magnitude 
of the measured cross section. 
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Figure 10: Detector-level distributions of reconstructed and calibrated Arj^ values for the mea¬ 
sured CG sample with a central LRG. The data are compared to predictions of (a) PYTHIA 
8 MBR, and (b) PYTHIA 8 4G simulations normalized to the integrated luminosity of the data 
sample. Gontributions for each of the generated processes are shown separately 


Table 2: The differential central pseudorapidity gap (DD dominated) cross section da/dAtj 
for Arj > 3, log^gMx > 1-1, and log^Q My > 1.1. The first and second errors correspond to 
statistical and systematic uncertainties, respectively. 


Ar] bin 

dace/dAtj (mb) 

3.0 < Arj < 4.5 

0.25 ± 0.003+°;°^ 

4.5 < At] < 6.0 

0.11 ±0.002+°;°3 

6.0 < At] < 7.5 

0.032 ± 0.001 ± 0.009 


8 Integrated cross sections 

The forward and central pseudorapidity gap samples are also used to measure the integrated 
cross sections in the kinematic regions given in Sections]^ andThe forward pseudorapidity- 
gap cross sections, ctVio-castor and Pcastor/ are measured in the region -5.5 < log^g < 
—2.5, for events without and with a GASTOR tag, corresponding to log^gMy < 0.5, and 
0.5 < log^^gMy < 1.1, respectively, while the central pseudorapidity-gap cross section, acc^ 
is measured for At] > 3, logjg Mx >1.1, and log^g My >1.1. 

Each cross section is evaluated by means of the formula 


a = 


^evt 

~AC' 


(5) 


where is the number of events in the kinematic regions given above, A is the acceptance, 
defined as the ratio of the number of events reconstructed to the number of events generated 
in that bin, taking into account the pileup correction, and C is the integrated luminosity. The 
acceptance is evaluated with the PYTHIA 8 MBR MG generator. 


Values of cr, 


I-CASTOR ^ 2.99 ± 0.02 (stat)+^ § (syst) mb, itcastor = 1.18 ± 0.02 (stat) ± 0.13 (syst) mb, 

^+0.13 


and acG = 0.58 ± 0.01 (stat)_Qj'j (syst) mb are obtained. Systematic uncertainties are evaluated 
as discussed in Section 11 As a consistency check of the analysis procedure, we measure the 
part of the total inelastic cross sections that is visible in the central GMS detector, correspond¬ 
ing to the region logjg Mx > 1.1 or logjg My > 1.1. A value of = 61.29 ± 0.07 (stat) mb 
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8 Integrated cross sections 
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Figure 11: The central pseudorapidity-gap cross section du /dArj (DD dominated) compared to 
MC predictions: (a) PYTHIA 8 MBR, PYTHIA 8 4C, and PYTHIA 6 Z2*, and (b) PHOJET, QGSJET- 
II 03, QGSJET-II 04, and EPOS. Error bars are dominated by systematic uncertainties, which are 
discussed in Section [TT| 


is found, in good agreement with the published CMS result Pinei(^ > 5 x 10“^) = 60.2 ± 
0.2 (stat) ± 1.1 (syst) ± 2.4 (lumi) mb, measured in a slightly different kinematic region (Mx or 
My > 16.7GeV) 130 1. According to PYTHIA 8 MBR, the phase space difference between the two 
measurements corresponds to > 5 x 10“^) = 0.5 mb. 


Table 3: Measured Gno-CASTOR/ gcastor/ and acc cross sections, compared to predictions of MC 
models. The first and the second errors in the data correspond to statistical and systematic 
uncertainties, respectively. 


Cross section 

Gno-CASTOR(nab) 

SD dominated 

gcastor (nab) 
DD dominated 

GcG (mb) 

DD dominated 

Data 

2.99 ± 

1.18 ±0.02 ±0.13 

0.58 ± 0.01+°;j^ 

PYTHIA 8 MBR 

3.05 

1.24 

0.54 

PYTHIA 8 4C 

3.31 

1.10 

0.40 

PYTHIA 6 Z2* 

3.86 

1.52 

0.78 

PHOJET 

3.06 

0.63 

0.32 

QGSJET-II 03 

2.63 

0.48 

0.22 

QGSJET-II 04 

1.70 

0.78 

0.37 

EPOS 

2.99 

0.85 

0.31 


Table 1^ lists the measured cross sections together with the absolute predictions of the MC sim¬ 
ulations. Based on the results presented thus far, the following conclusions can be drawn 
about the models: PHOJET, QGSJET-II 03, QGSJET-II 04, and EPOS predict too few DD events, 
which dominate the measured Pcastor and Gcg cross sections (Figs. and [TT]?); among 
these four models only QGSJET-II 03 satisfactorily predicts the Gno-CASTOR cross section (Fig.|^); 
PYTHIA 8 4C, PYTHIA 6 Z2*, PHOJET, and EPOS do not predict correctly the dependence for 
the SD process, which dominates the measured forward pseudorapidity gap cross section for 
log^Q My < 0.5 (Fig.[^a,b); and PYTHIA 8 MBR describes the data within uncertainties in all the 
measured regions. 
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9 The SD and DD cross sections 

The c^no-CASTOR cross section discussed above is dominated by SD events (Fig. [S]?), whereas 
the cTcASTOR arid acQ cross sections are mainly due to DD events (Figs. and |l0|. As the 
contribution from ND and other diffractive processes to these cross sections is small, we use 
the event decomposition as defined in the pythia 8 MBR simulation with e = 0.08 to correct 
for them and extract the SD and DD cross sections. 

The dominant background in the cyio-cASTOR cross section originates from DD events; the CD 
contribution is minimal, while the ND contribution is negligible (Fig.|^). The DD contribution 
is well understood via the CASTOR-tag events (Fig. |^), and has an uncertainty of ~10-20% 
due to the ND contamination. Since the DD events contribute about 20% to the no-CASTOR- 
tag sample, the uncertainty in the SD cross section due to the subtraction of the DD component 
amounts to only a few percent. The visible part of the total SD cross section, corresponding 
to —5.5 < log^g < —2.5, is found to be = 4.06 ± 0.04(stat)lQ;g3 (syst)mb. The result 

accounts for both pp ^ Xp and pp —> pY. 

The dominant background to the (Tcastor ^rid ctcg cross sections originates from ND events. 
The CD and SD contributions are negligible in the CASTOR-tag sample, while SD events con¬ 
tribute minimally to the central-gap sample (Figs. and [lO^). The total DD cross sections 
integrated over the regions 

• —5.5 < logjQ < —2.5 and 0.5 < log^g My < 1.1, and 

• At] > 3, logjQ Mx > 1.1 and log^g My > 1.1, 

are c^castor “ ^ (stat) ± 0.12 (syst) mb, and ~ ^ i 3 

respectively. 

To provide the DD cross section in the widest kinematic region spanned by the data, we also 
evaluate the visible DD cross section, defined as = 2 g^astor where the 

factor of 2 assumes the same dependence of the DD cross section on Mx and My (boxed regions 
in Fig.|^). This leads to = 2.69 ± 0.04 (stat)]'^Q 3 Q (syst), in the kinematic region delimited 

by the solid and dashed (red) lines in Figs. |^,c. This result is used below to extrapolate the DD 
cross section to the region Arj > 3 (to the left of the dashed (blue) line in Fig.|^). 

9.1 Extrapolation of the visible SD and DD cross sections 

The measurements based on the central CMS detector are insensitive to the low-mass part 
of diffractive dissociation. Therefore, in order to compare the measured cross section 

with results of other experiments and theoretical models that present integrated cross sections 
for ^ < 0.05, an extrapolation from —5.5 < log^g ^x < —2.5 (^x = M^/s) to ? < 0.05 is 
required. Similarly, the cr^^vis ^j-oss section must be extrapolated to Arj > 3 (Fig. The ex¬ 
trapolation factors, calculated by using each of the MC simulations introduced in Section]^ are 
presented in Appendix A. Not all of the simulations are able to describe the measured cross 
sections (see Section |^, nor do they include realistic hadronization models (see Appendix 
B). Following the discussion in Section]^ and Appendix B, the extrapolation factors are de¬ 
termined with PYTHIA 8 MBR (with e = 0.08), which describes well all aspects of our data. 
The multiplicative factor needed to extrapolate the measured SD cross section to ^ < 0.05 
is = 2.18^4yJ°, and that for the extrapolation of the DD cross section to Arj >3 is 

/mbr — (Tables 1^ and in Appendix A). The extrapolation uncertainties are esti¬ 

mated by changing the parameters a' and e of the Pomeron trajectory from their nominal 
values (a' = 0.25 GeV^ , e = 0.08) to those presented in Tables and (one parameter 
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9 The SD and DD cross sections 


changed at a time), and adding in quadrature the corresponding deviations with respect to 
the central result, separately for the positive and negative deviations. The extrapolated SD 
and DD cross sections thus obtained are = 8.84 ± 0.08 (stat)+J gg (syst)^Q 37 (extrap) mb and 
= 5.17 ± 0.08 (stat)lQ 57 (syst))')Q 33 (extrap) mb, respectively. 




Figure 12: Diffractive cross sections as a function of collision energy measured in pp and pp 
collisions l3Tti36l compared to pythia 8 MBR (e = 0.08,0.104) and other model predic¬ 
tions ll37H3^ : (a) total SD cross section for ^ < 0.05, and (b) total DD cross section for Arj > 3. 
The irmer (outer) error bars of the CMS data points correspond to the statistical and systematic 
(and the additional extrapolation) uncertainties added in quadrature. 


Figure [l^ presents the extrapolated SD cross section compared to the ALICE result and a 
compilation of lower center-of-mass energy measurements ll3TM35ll . The data are also compared 
to the PYTHIA 8 MBR simulation, as well as to the GLM Il37 l and KP Il38l l39l models. The CMS 
result is consistent with a SD cross section weakly rising with energy. 


Figure!!^ shows the extrapolated DD cross section compared to the ALICE results lH, those by 
CDF atVs = 630 GeV and 1.8 TeV BMl, as well as the pythia 8 MBR, GLM EZi, and KP ||38ll^ 
models. The GMS measurement at = 7 TeV is in agreement with the ALIGE measurement at 




















9.2 Summary of results 
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Table 4: Measured diffractive cross sections in regions of Mx (or ^x = 
= My/s), and Arj {At] = — log^, where ^ = M|'My/(swjp) for DD). 
for the cross section extraction is indicated as LRG for calculations involving 

- M|;/s), My (or 

The method used 
: all events selected 

in the LRG sampl 

es, and as MBR for calculations that involve background subtraction or ex- 

trapolation based 

on the prediction of the PYTHIA 8 MBR simulation. The first and the second 

errors in the data correspond to the statistical and systematic uncertainties, respectively For 
and (7^^, the third errors correspond to the extrapolation uncertainties. 

Gross section 

Mx or ^x range 

My or ^Y range 

Ai] range 

Result (mb) 

LRG 





Gio-CASTOR 

-5.5 < log^x < -2.5 

logio < d.5 

— 

2.99 ± Omtoil 

Gastor 

-5.5 < log^x < -2.5 

0.5 < log^Q My < 1.1 

— 

1.18 ±0.02 ±0.13 

Gg 

log^^gMx >1.1 

logioMy >1.1 

At] > 3 

0.58 ± 0.0i;°;|3 

MBR 





g-SDvis ~ 

-5.5 < log^x < -2.5 

Mx = nip 

My = nip 

—5.5 < log^y < —2.5 

— 

4.06 ± 0.04^°;^^ 

„-DDvis 

‘^CASTOR 

-5.5 < log^x < -2.5 

0.5 < log^Q My < 1.1 

— 

1.06 ±0.02 ±0.12 

,,-DDvis 

Gcg 

log;^oMx >1.1 

log^gMy >1.1 

At] > 3 

0.56 ± 0.0i;|];|^ 

-5.5 < log^x < -2.5 
0.5 < logio Mx < 1.1 
log^^gMx >1.1 

0.5 < log^Q My < 1.1 
—5.5 < log^y < —2.5 

log^oMy >1.1 

At] > 3 

2.69 ± 0.04;0;2^ 

MBR 





^SD 

< 0.05 

Mx = nip 

My = nip 
< 0.05 

— 

8.84 ± 0.08;j;49 

^DD 

— 

— 

At] > 3 

5.17±0.08l°f 


the same energy. Note, however, that the ALICE result is obtained from the NSD (non-single- 
diffractive = DD+ND) data, while for the CMS measurement the ND background has been 
subtracted. Here as well, the data are consistent with a weakly rising cross section with energy, 
as predicted by the models. 

9.2 Summary of results 

Table presents the summary of the cross section measurements illustrated in the previous 
sections, together with the kinematic region covered by each measurement. The method used 
for the cross section extraction (LRG or MBR) is given as well. The cy,o-CASTOR/ Gzastor/ arid 
acQ cross sections are measured from all the events passing the LRG selection (Section |^. The 
^SDvis, and cross sections are extracted from the latter ones by subtracting the 

background contribution from other processes as predicted by the pythia 8 MBR simulation 
(Section]^. In addition, (jr’Dvis calculated from the combination of c^castor ^^d 
nally, the and cross sections (Section |9.l[ ) are calculated by extrapolating c^rivis 

crDDvis |.Q region of lower diffractive masses using the mass dependence of the cross section 
predicted by PYTHIA 8 MBR. 
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10 Pseudorapidity gap cross section 


10 Pseudorapidity gap cross section 

This section presents the results of an alternative approach to the study of diffractive events, in 
which the data are analyzed in terms of the widest pseudorapidity gap adjacent to the edge of 
the detector BOl . In each event, particles are first ordered in t], and the largest pseudorapidity 
gap, , is determined as ls,rf = max(|?/n\m — |?/max — where = ±4.7 are the 
detector edges in r], and r]max is the highest (lowest) rj of the PF objects in the event (see 

Fig.§. 
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Figure 13: Uncorrected Arj^ distribution compared to various detector-level MC predictions. 


The analysis is based on a minimum bias data sample, selected as described in Section with 
negligible pileup (0.006), and corresponding to an integrated luminosity of 20.3 The un¬ 
corrected distribution of the pseudorapidity gap size is shown in Fig. along with the pre¬ 
dictions of various MC models. A wider bin width is used at low Arj^ to account for the lower 
spatial resolution in the forward region. 


10.1 Corrections for experimental effects 

Interactions of the beam protons with the residual gas particles in the beam pipe or inside the 
detector region affect the pseudorapidity gap distribution in data. The overall beam-induced 
background, integrated over the full measurement region, is about 0.7%. After the subtrac¬ 
tion of this background, the differential cross section dcr/ dAr]^ is determined according to the 
formula 


da _ 

dAl/F Ui: (Al/F)bin' 

where N®'"* is the number of events in the bin, corrected for migration effects, U the trigger 
efficiency, C the integrated luminosity, and (A?^^)bin the bin width. 

The trigger efficiency is obtained from a comparison with zero-bias data where no trigger re¬ 
quirements were applied. In order to have a satisfactory correlation between the generated and 
reconstructed values of Arj^, and hence a reliable correction for bin-migration effects, the cross 















10.2 Corrected results 


19 


section is evaluated for events with at least one stable final-state particle of transverse momen¬ 
tum pj > 200 MeV in the region of \t]\ <4.7. The migration corrections are evaluated with the 
iterative Bayesian unfolding technique Il28]l . as implemented in the RooUnfold package Il29l . 
with a single iteration. The response matrix is obtained with PYTHIA 8 MBR (e = 0.08). The 
cross section is measured only for At]^ < 8.4, so as to avoid regions where the trigger ineffi¬ 
ciency and the unfolding uncertainty are large. 


10.2 Corrected results 

Figure [^shows the unfolded and fully corrected differential cross section da/dArj^ for events 
with at least one particle with pj > 200 MeV in the region of |j/| < 4.7. As the statistical 


uncertainty is negligible, only the systematic uncertainty, discussed in Section 11 is shown. 
The predictions from PYTHIA 8 MBR (e = 0.08 and 0.104), PYTHIA 8 tune 4C, and PYTHIA 6 
tune Z2* are also given. The MC predictions show that in the pseudorapidity range covered 
by the measurement, \t]\ < 4.7, a large fraction of nondiffractive events can be suppressed by 
means of the Arj^ > 3 requirement. 

The present results are consistent with those from the ATLAS collaboration |^, as shown in 


Fig. 15 The ATLAS measurement uses all stable final-state particles with pj > 200 MeV over 
the region \f]\ < 4.9. According to the PYTHIA 8 MBR simulation, the difference in the t] cover¬ 
age between the two experiments causes changes in the Arj^ distribution of up to 5%. 


11 Systematic uncertainties 

Systematic uncertainties are obtained by varying the selection criteria and modifying the anal¬ 
ysis. The following sources of systematic uncertainties are taken into account for the results 
presented in Sections [6}jT0{ 

• HF energy scale: varied in the MC simulations by ±10%, to reflect the energy scale 
uncertainty estimated for data. 

• PF energy thresholds: raised by 10%, based on dedicated studies of the detector 
noise. 

• Modeling of the diffractive interaction and the hadronization process: the hadroni- 
zation parameters in the nominal PYTHIA 8 MBR MC sample are tuned to describe 
the multiplicity and pj spectra of diffractive systems in pp and pp collisions at 
^/s < 1800 GeV |4Til . The corresponding uncertainty is estimated by taking the 
difference (see Appendix B) between the results obtained with PYTHIA 8 MBR and 
those obtained with PYTHIA 8 4C (Sections | 6 ]-[^, or PYTHIA 8 4C and PYTHIA 6-Z2* 
(Section [10|. 

• The uncertainty in the integrated luminosity measurement is ±4% (421 113 . 

In addition, the following checks are carried out for the results shown in Sections | 6 ]-^ 

• CASTOR energy scale ||44|| : changed in the simulation by ±20%, to reflect the esti¬ 
mated energy scale uncertainty for the data. 

• CASTOR energy threshold in each sector fl^ : changed from the nominal Aa to 3.5a 
and 5a, where a is the pedestal width. 

• CASTOR alignment uncertainty ||45]| : the simulated CASTOR position in the plane 
transverse to the beamline is varied within the limits allowed by the condition that 
the MC description of the energy flow in CASTOR remains satisfactory This corre- 
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11 Systematic uncertainties 


20.3 (7 TeV) 20.3 (7 TeV) 




Figure 14: Differential cross section dcr/for stable particles with pj > 200 MeV in the 
region \r]\ < 4.7 compared to the corresponding predictions of (a) PYTHIA 8 MBR (e = 0.08), (b) 
PYTHIA 8 MBR (e = 0.104), (c) PYTHIA 8 tune 4C, and (d) PYTHIA 6 tune Z2*. The band around 
the data points represents the total systematic uncertainty, which is discussed in Section 
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sponds to about 10 mm and 4 mm, for the left and right CASTOR sides, respectively. 

• Trigger efficiency uncertainty: estimated from a comparison of efficiency curves be¬ 
tween data (measured by using a control sample for which no trigger requirements 
were applied) and MC. 

• Background subtraction: backgrounds from DD and ND events in the SD sample, 
and from ND and SD in the DD sample are estimated with PYTHIA 8 MBR (Figs. 
and[T0|. The corresponding uncertainty is estimated by varying their relative contri¬ 
butions by ±10% (average normalization uncertainty of the model). The contribu¬ 
tion from CD events in the SD and DD samples is negligible. 
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20.3 nb-' (7 TeV) 
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Figure 15: Differential cross section dcr/for stable particles with pj > 200 MeV in the 
region \ri\ < 4.7 compared to the ATLAS result [JJ: distributions (top) and ratio to the CMS 
measurement (bottom). The band represents the total systematic uncertainty in the CMS mea¬ 
surement, while the uncertainty in the ATLAS measurement is shown by the error bars. The 
stable-particle level definitions of the two measurements are not exactly identical: CMS mea¬ 
sures the forward pseudorapidity gap size starting from p = ±4.7, whereas the ATLAS limit is 
rj = ±4.9. 


The total systematic uncertainty is obtained by summing all individual uncertainties in quadra¬ 
ture, separately for the positive and negative deviations from the nominal cross section values, 
leading to a total systematic uncertainty of up to 25%. Tablepresents the summary of the 
systematic uncertainties for the measurement of the cross sections. 

The systematic uncertainties are significantly larger than the statistical ones, and the dominant 
sources are the HF energy scale and the modeling of diffraction and hadronization. For the 
^CASTOR section, also the uncertainty related to the CASTOR alignment is significant. 


12 Summary 

Measurements of diffractive dissociation cross sections in pp collisions at = 7 TeV have 
been presented in kinematic regions defined by the masses Mx and My of the two final-state 
hadronic systems separated by the largest rapidity gap in the event. Differential cross sections 
are measured as a function of ^x = ^x/® region —5.5 < log^o < -2.5, for log^o < 
0.5, dominated by single dissociation (SD), and 0.5 < log^gMy < 1.1, dominated by dou¬ 
ble dissociation (DD). The discrimination between the above two My regions is performed by 
means of the CASTOR forward calorimeter. The cross sections integrated over these regions are 
Vio-CASTOR = 2.99 ± 0.02 (stat)lQ 29 cVastor = 1-18 ± 0.02 (stat) ± 0.13 (syst) mb, 

respectively. 

The inclusive pp cross section is also measured as a function of the width of the central pseu¬ 
dorapidity gap, Arj, for Arj > 3, log^Q Mx > 1-1, and log^Q My > 1.1 (dominated by DD contri¬ 
butions). The corresponding integrated cross section is gcg = 0.58 ± 0.01 (stat)^Qj^ (syst) mb. 
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12 Summary 


Table 5: Systematic uncertainties for the measurement of the ^castor' cross 

sections; individual contributions, as well as total systematic and statistical uncertainties are 
shown. 


Source 

Uncertainty (%) 

^SDvis 

^DDvis 

^CASTOR 

,T-DDvis 

•^CG 

HE energy scale 

10 

1.6 

23 

PE thresholds 

0.8 

0.4 

6.9 

Diff. and had. modeling 

10 

4.3 

0.4 

Luminosity 

4 

4 

4 

CASTOR energy scale 

0.5 

0.9 

0 

CASTOR threshold 

0.9 

2.8 

0 

CASTOR alignment 

2.6 

8.3 

0 

Trigger 

0.6 

0.6 

0.7 

Background sub. 

4.3 

0.4 

1.3 

Total systematic 

16 

11 

25 

Statistical 

0.9 

1.8 

1.3 


The SD and DD cross sections in the above three regions, extracted by means of the PYTHIA 
8 MBR simulation, which provides a good description of the data, are = 4.06 ± 0.04 (stat)+Q ^3 (syst) mb 

(accounting for both pp —> Xp and pp —> pY), = 1.06 ± 0.02 (stat) ± 0.12 (syst) mb, and 

= 0.56 ± 0.01 (stat)+Q ([3 (syst) mb, respectively. 

Extrapolations of the SD and DD cross sections to the regions ^ < 0.05 and A// > 3, performed 
with PYTHIA 8 MBR, yield ( 7 ®^ = 8.84 ± 0.08 (stat)+{ 33 (syst)lQ 37 (extrap) mb and = 5.17 ± 

0.08 (stat)lQ 57 (syst)('( q 53 (extrap)mb, respectively. 

In addition, the inclusive differential cross section d( 7 /dA//^ for events with a pseudorapidity 
gap adjacent to the edge of the detector is measured over 8.4 units of pseudorapidity. 

These measurements are compared to results from other experiments as well as to phenomeno¬ 
logical predictions. The data are consistent with the SD and DD cross sections weakly ris¬ 
ing with energy, and provide new experimental constraints on the modeling of diffraction in 
hadronic interactions. 
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A SD/DD extrapolation factors 

Figure [l^ shows the = M\/s dependence of the SD cross section for the pythia 8 4C, 
PYTHIA 6 72*, PHOJET, QGSJET-II 03, QGSJET-II 04, and EPOS simulations, compared to the nom¬ 
inal PYTHIA 8 MBR simulation used in this analysis for two regions of ^x, —5.5 < log^Q < 
—2.5 (dashed yellow) and ^x < 0.05 (solid khaki). In addition, the pythia 8 MBR simulations 
with values of a' and e changed to a' = 0.125 GeV^^, e = 0.104, and e = 0.07 (one parameter 
changed at a time) are also included to provide a scale for their effect on the cross sections. Ex¬ 
trapolation factors, defined as the ratios of e‘®’^(^x < 0.05) to g®''"*(— 5.5 < logio^x < -2.5), 
are presented for each of the above ten MC models in Table For each model, two ratios are 
evaluated, one in which both cross sections (numerator and denominator of the extrapolation 
factor) are calculated by using the same generator (/®), and another where the prediction of 



for PYTHIA 8 4C, PYTHIA 6 72*, PHOJET, QGSJET-II 03, QGSJET-II 04, EPOS MC, and PYTHIA 
8 MBR with the parameters of the Pomeron trajectory changed from the nominal values (a' = 
0.25 GeV^^, £ = 0.08) to a' = 0.125 GeV^^, e = 0.104, and e = 0.07 (one parameter changed at a 
time). The nominal PYTHIA 8 MBR simulation is presented in each plot for the two regions of 
^x, —5.5 < log^g ^x < —2.5 (dashed yellow) and ^x < 0.05 (solid khaki), used to extrapolate 
the measured SD cross section (from the dashed (yellow) to the solid (khaki) regions). 
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A SD/DD extrapolation factors 


the 5.5 < log^Q^x < —2.5) is calculated W using the nominal pythia 8 MBR with 

£ = 0.08 (/]^r). The numbers in brackets in Table 0 show the relative change of the extrapo¬ 
lation factors with respect to the nominal PYTHIA 8 MBR simulation; the one related to is 
sensitive to the difference in the shape of the (mass) distribution, while the one related to 
/]^R sensitive to the the normalization of the SD cross section. 

Table ^ shows the extrapolation factors for the DD case, defined as the ratios of a^^{Ar] > 
3) to again for two cases: one when both cross sections are calculated with the same 

MC generator (/^^), and the other when the predicted ( 7 ^d,vis ^j-qss section is from PYTHIA 
8 MBR with e = 0.08 (/mbr)- The relative change of the extrapolation factors with respect to 
the nominal PYTHIA 8 MBR simulation is shown in brackets; the one related to accounts 
for the difference in the shape of the Mx and My dependence, while the one related to /mbr 
is sensitive to the difference in the mass dependence and the normalization of the DD cross 
section. 

Table 6: Extrapolation factors= C7f°(^ < 0.05)and= C7f°(^ < 0.05) 
from the visible to total SD < 0.05) cross section for each MC model considered (f = 1 — 10). 
The relative change with respect to the value obtained by PYTHIA 8 MBR with e = 0.08 is shown 
in parenthesis. 


i 

MC model 

jrSD 

fSD 

iMBR 

1 

PYTHIA 8 MBR (£ = 0.08) 

2.18(1.00) 

2.18(1.00) 

2 

PYTHIA 8 4C 

2.32(1.06) 

2.51(1.15) 

3 

PYTHIA 6 Z2* 

2.29(1.06) 

2.89(1.34) 

4 

PHOJET 

2.06 (0.95) 

2.18(1.00) 

5 

QGSJET-II 03 

2.72(1.25) 

3.19(1.46) 

6 

QGSJET-II 04 

3.62(1.66) 

2.30 (1.06) 

7 

EPOS 

3.44(1.58) 

2.15(0.99) 

8 

PYTHIA 8 MBR (a' = 0.125) 

2.27(1.04) 

2.34(1.07) 

9 

PYTHIA 8 MBR (£ = 0.104) 

2.23(1.03) 

2.42 (1.11) 

10 

PYTHIA 8 MBR (£ = 0.07) 

2.16(0.99) 

2.09 (0.96) 


Table 7: Extrapolation factors = aP^{At] > 3)/i7'^°''''® and /^br = > 3)/(7^g^^® 

from the visible to total DD {Atj > 3) cross section for each MC model considered (i = 1 — 10). 
The relative change with respect to the value obtained by PYTHIA 8 MBR with e = 0.08 is shown 
in parenthesis. 


i 

MC model 

yDD 

rDD 

Imbr 

1 

PYTHIA 8 MBR (e = 0.08) 

1.92 (1.00) 

1.92 (1.00) 

2 

PYTHIA 8 4C 

2.52(1.32) 

1.86 (0.97) 

3 

PYTHIA 6 Z2* 

2.39(1.25) 

2.15 (1.13) 

4 

PHOJET 

1.80 (0.94) 

0.60 (0.31) 

5 

QGSJET-II 03 

— 

— 

6 

QGSJET-II 04 

2.04(1.07) 

0.94 (0.49) 

7 

EPOS 

4.73 (2.47) 

1.93(1.01) 

8 

PYTHIA 8 MBR (a' = 0.125) 

1.97(1.03) 

2.32(1.21) 

9 

PYTHIA 8 MBR (£ = 0.104) 

2.00(1.04) 

2.37(1.24) 

10 

PYTHIA 8 MBR (e = 0.07) 

1.88 (0.98) 

1.73(0.90) 
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B MC hadronization models 

In this section, the hadronization models used to generate particle spectra in the simulations 
introduced in Section]^ are compared to a reference model 1^1461 based on data. The model 
correctly describes the charged-particle multiplicity and pj spectra of diffractive and inclusive 
proton-(anti)proton data for \/s < 1800 GeV by assuming that the Pomeron-proton collision 
produces a system of mass Mx that hadronizes as if it had been produced in a nondiffractive 
proton-proton collision at ^/s = Mx- Figures]^ andshow the charged-particle multiplicity 
distributions and pj spectra for the SD process for three ranges of Mx'- 5.6 < Mx < 10 GeV, 
32 < Mx < 56 GeV, and 178 < Mx < 316 GeV, for the pythia 8 MBR, pythia 8 4G, PYTHIA 
6 Z2*, PHOJET, QGSJET-II 03, QGSIET-II 04, and EPOS simulations, compared to the reference 
model quoted above. The following conclusions can be drawn: PYTHIA 6 Z2*, QGSJET-ll 03, 
QGSJET-II 04, and EPOS predict smaller multiplicities and harder pj spectra than the model of 
Ref. ||4T1 l46ll : pythia 8 4G agrees with the mean values of the multiplicity distributions, but 
predicts narrower widths, and harder pj spectra; PHOJET predicts multiplicity distributions 
consistent with the reference model with harder pj spectra; and PYTHIA 8 MBR agrees with 
the reference model in both multiplicity and px spectra. The latter is an expected result, as 
the hadronization parameters of the pythia 8 MBR simulation have been tuned to follow the 
reference model. 
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B MC hadronization models 
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Figure 17: Charged-particle multiplicity (Nch) distributions (area-normalized) in the PYTHIA 
8 MBR, PYTHIA 8 4C, PYTHIA 6 Z2*, PHOJET, QGSJET-II 03, QGSJET-II 04, and EPOS MC simula¬ 
tions (rows) in three bins of Mx (columns) in SD collisions, compared to a reference hadroni¬ 
zation model (dashed line), which describes the available data at \/s < 1800 GeV ll4Tll46l . 
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Figure 18: Transverse-momentum (pj) distributions (area-normalized) in the PYTHIA 8 MBR, 
PYTHIA 8 4C, PYTHIA 6 Z2*, PHOJET, QGSJET-II 03, QGSJET-II 04, EPOS MC simulations (rows) 
in three bins of Mx (columns) in SD collisions, compared to a reference hadronization model 
(dashed line), which describes the available data at < 1800 GeV ETlIi^ . 
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